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Developing a theory of high-temperature superconductivity in copper oxides is one 
of the outstanding problems in physics. Twenty-flve years after its discovery [1], no 
consensus on the microscopic theory has been reached despite tremendous theoretical 
and experimental efforts. Attempts to understand this problem are hindered by the 
subtle interplay among a few mechanisms and the presence of several nearly degener- 
ate and competing phases in these systems. Here we provide unified parameter-free 
explanation of the observed oxygen-isotope effects on the critical temperature, the 
magnetic-field penetration depth, and on the normal-state pseudogap for underdoped 
cuprate superconductors within the framework of the bipolaron theory compatible 
with the strong Coulomb and Prohlich interactions, and with many other indepen- 
dent observations in these highly polarizable doped insulators. Remarkably, we also 
quantitatively explain the measured critical temperature and the magnitude of the 
magnetic-field penetration depth. The present work thus represents an ultimate proof 
of the bipolaron theory of high-temperature superconductivity, which takes into ac- 
count essential Coulomb and electron-phonon interactions. 



In 1911 curiosity concerning the electrical properties 
of metals at low temperatures led the Dutch physicists, 
Kamerling Onnes and his assistant G. Hoist to discover 
superconductivity at 4.2 K in mercury This discov- 
ery was one of the most important experimental find- 
ings in low temperature physics. On the long way to- 
wards a microscopic understanding of superconductivity, 
the observation of an isotope effect on the critical tem- 
perature, Tc in 1950 [3, 4] gave an important clue to 
the microscopic mechanism of superconductivity. The 
presence of an isotope effect thus implies that supercon- 
ductivity is not of purely electronic origin. In the same 
year Frohlich [5] pointed out that the electron-phonon 
interaction gave rise to an attractive interaction between 
electrons, which might be responsible for superconduc- 
tivity. Frohlich's theory played a decisive role in estab- 
lishing the correct mechanism. Finally, in 1957, Bardeen, 
Cooper and Schrieffer \7j (BCS) developed the BCS the- 
ory that was the first successful microscopic theory of 
superconductivity. The BCS theory implies an isotope- 
mass dependence of with an isotope-effect exponent 
a ~ — din Tc /din A/ = 1/2, in excellent agreement with 
the reported isotope exponents in simple metallic super- 
conductors like Hg, Sn and Pb. 

The doping dependent oxygen-isotope effect (OIE) 
on the critical temperature Tc, oP — — din Tc/d In Mo 
(where Mq is the oxygen-isotope mass) [8] and the sub- 
stantial OIE on the in-plane su per carrier mass fn*^, 
a^, — dm*^/d\nMo (ref. 0-14|), provide direct 
evidence for a significant electron-phonon interaction 
(EPI) also in high-temperature cuprate superconduc- 
tors. High resolution angle-resolved photoemission spec- 
troscopy (ARPES) [l3| provides further evidence for 



the strong EPI apparently with c-axis-polarised optical 
phonons. These results along with optical jl6| . neutron 



scattering 17, and tunneling data 13-21 
biguously show that lattice vibrations play a significant 
but unconventional role in high-temperature supercon- 
ductivity. The interpretation of the optical spectra of 
high-Tc materials as due to multi-polaron absorption 22 1 
strengthens the view [2^ that the Frohlich EPI is im- 
portant in those structures. Operating together with a 
shorter-range deformation potential and molecular-type 
(e.g., Jahn-Teller [22|) EPIs, the Frohhch EPI can read- 
ily overcome the Coulomb repulsion at a short distance 
of about the lattice constant for electrons to form real- 



space intersite bipolarons (see 25| and Supplementary 



Information 26[) 



Despite all these remarkable and well-done experi- 
ments that lead to the consistent conclusion about the 
important role of EPI in high-temperature superconduc- 
tors, there is no consensus on the microscopic origin of 
the observed unconventional isotope effects on the in- 
plane magnetic-field penetration depth and the normal- 
state pseudogap. The doping dependent a*^ has been 
explained as due to the doping independent oxygen- 
isotope effect on the in-plane carrier concentration n, 
that is, a° = -dlnn/dlnMo = 0.146 (ref. ^). This in- 
terpretation contradicts other independent experiments 



terp r 

[111 Il2l Il4j which consistently show that the carrier con- 
centrations of the two oxygen-isotope samples are the 
same within 0.0004 per Cu. This is also in sharp con- 
trast to the observed very large oxygen-isotope effect on 
the low-temperature magnetic-field penetration depth, 
Kb oc {m**Jnf/^, in both Lai.94Sro.o6Cu04 (ref. [123) 
and Yo.55Pro.45Ba2Cu307_y (ref. 28] ), which would lead 
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FIG. 1: (Color online) Low energy excitations in cuprate su- 
perconductors. 



to ~ 2 if one would assume that the supercarrier mass 
is independent of the oxygen mass. Another model based 
on the pair-breaking effects due to impurities, disorder, 
and/or pseudogap can also explain the observed oxygen- 
isotope effects on the penetration depth and the critical 
temperature in deeply underdoped samples [l^ . But this 
model cannot consistently explain the negligibly small 
but the large OIE on the penetration depth in optimally 
doped samples jlS,] . 

Alternatively, the bipolaron theory of superconductiv- 
ity [2^ can naturally account for the substantial a^. and 
large in deeply underdoped cuprates js^]. There is a 
qualitative difference between ordinary metals and po- 
laronic conductors. The renormalized effective mass of 
electrons is independent of the ion mass M in ordinary 
metals (where the Migdal adiabatic approximation is be- 
lieved to be vahd), because the EPI coupling constant A 
does not depend on the isotope mass (see, for instance 
ref. [2^). However, when electrons form polarons (new 
quasiparticles dressed by lattice distortions), their effec- 
tive mass m* depends on M [sot . 

Although the bipolaron theory can qualitatively ex- 
plain both a^j, and a'-' in deeply underdoped samples, 
some important issues have not been well addressed by 
the theory. The first issue is why a^* is not equal to a'-' 
even in deeply underdoped cuprates. The second issue 
is why a*^ is much smaller than a^. for slightly under- 
doped samples. The third issue is why there is a giant 
negative oxygen-isotope effect on the pesudogap forma- 
tion temperature T* in HoBa2Cu408 (ref. [31i)- Here 
we provide parameter-free explanations to the observed 
oxygen-isotope effects on the critical temperature, the in- 
plane supercarrier mass, and on the normal-state pseu- 
dogap in HoBa2Cu408 within the framework of the bipo- 
laron theory. The present work thus represents an ulti- 
mate proof of the bipolaron theory of high-temperature 
superconductivity. 



We adopt the bipolaronic low-energy excitation struc- 
ture of cuprates, Fig. 1, derived from the microscopic 
Hamiltonian with the strong Coulomb and Frohlich in- 
teractions [25l.[26|. The critical temperature of quasi- two- 
dimensional (2D) bipolarons, which are hard-core bosons, 
depends on their density ni,(Tc) at the critical tempera- 
ture and the in-plane bipolaron mass Jn*^: 



Tr oc 



(1) 



"ab 



The bipolaron density slightly depends on temperature 
due to bipolaron depletion into unbound single polarons, 
nb{T) — [x ~ 7ip(T)]/2, with the density Up given by 



-A/2kBT 



(2) 



where W is the polaron half-bandwidth, x is the in-plane 
doping level, and A is the bipolaron binding energy. 
Fig. 1. This expression is obtained by integrating the 
Fermi-Dirac distribution function with a constant (2D) 
density of states in the polaron band, N{E) = 1/2W 
and assuming that the polaron half-bandwidth is large 
enough, W > /csTc, 



np{T) 



2W 



dE 



2N{E) 



l + exp[(A/2 + £;)/fcBr] 



(3) 



Here we take into account that the chemical poten- 
tial is zero in the superconducting state at and below 
Tc if all energies are taking with respect to the bipo- 
laron ground state. The polaron in-plane effective mass. 



m*^ oc exp(ylvM) {A is a constant), the polaron inverse 
bandwidth, 1/W, and the inter-site bipolaron mass m*^ 
have the same isotope exponent dlnTO*^/lnM = am* = 
(1/2) ln(TO*j/TO) [30|, where m is the band mass in a rigid 
lattice. The isotope effect on the pseudogap is given by: 
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(4) 



(5) 



The above expressions are obtained by taking into ac- 
count that the pseudogap in Fig. 1 is given by A/ 2 — 
Jp/2-(W-t/2) = Jp/2-3W/A, where Jp is the phonon- 
induced intersite attraction, which is independent of the 
ionic mass [2^ , and t « W/ 2 in the intermediate coupling 
regime Then using Eqs. (HI [2|) and neglecting the 

terms on the order of kBTc/W <C 1, one readily obtains 
the ratio: 



= 1 - 



1 



2p(T,)][l+exp(A/2fcsr,)]' 



(6) 



Eq. [6] can naturally explain why aj^* is always larger 
than [ll|, ■ It is worth noting that Eq. 11] is valid 
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FIG. 2: (Color online) The relaxation rate F of crystal-field 
excitations for the ^^0 and samples of slightly under- 
doped HoBa2Cu408. The data are taken from ref. [31(. The 
solid red lines represent the best fits of Eq. [8]to the data be- 
low T*. The best fits yield A/its = 94.6±8.8 K for the ^''O 
sample and A/ks = 286±19 K for the sample. 



only if din A/dln M is small. When c?lnA/dlnM is 
large, we need to use Eqs. [Hand [2] to directly calculate 
the Tc and np{Tc) changes upon the isotope exchange, 
that is, 



5T, 

T, 



(7) 



Using the above equations we can quantitatively ex- 
plain the oxygen-isotope effects on the pesudogap and the 
critical temperature in slightly underdoped HoBa2Cu408 
(ref. [31]). The oxygen isotope effect on the relaxation 
rate of crystal-field excitations in this compound was in- 
vestigated by means of inelastic neutron scattering I^H . 
The relaxation rate, which is related to the free-carrier 
spin density, clearly shows a large oxygen-isotope effect 
(see Fig. 2). For the "'^^O sample there is evidence for the 
opening of an electronic gap in the normal state at T* 
~ 170 K while for the sample T* is shifted to about 
220 K. In contrast, the is shifted from 79.0 to 78.5 K 
upon replacing ^^O with ^^O (the concentration is 
about 75%) [3l|. 

In a normal state with no pseudogap, the relaxation 



rate F„(r) is proportional to \Je.xN{E,p)'\^T ^ where J^x 
is the exchange integral between the 4/ electrons of the 
Ho'^"'" ions and the charge carriers and N(EpY& the elec- 
tronic density of states at the Fermi energy [3l| . Within 
the polaron/bipolaron framework, N{Ep) = 1/2W oc 
m*^ cx oxygen-isotope effect on N{Ep) is 

the same as the oxygen isotope effect on m*^. For a 
slightly underdoped YBa2Cu3 07_j^ film, Sm'^l/m'^l was 
found to be 5.5% upon replacing ^'^O with ^^O (the ^*^0 
concentration is about 95%) [14]. If we assume that the 
oxygen-isotope effect on m*^ for HoBa2Cu408 is similar 
to that for the slightly underdoped YBa2Cu307_y film, 
we expect that SmU/mH = 4.3% in HoBa2Cu408. 

In the superconducting state or below T*, the relax- 
ation rate is suppressed due to openning of the gap. Then 
the relaxation rate is given by ^] : 



F cx Texp(-A/2fcBr). 



(8) 



The solid lines in Fig. 2 represent the best fits of 
Eq.[S]to the data below T*. The best fits yield A/ks 
94.6±8.8 K for the ^''O sample and A/fcs = 286±19 K 
for the ^^O sample. Therefore, there is a giant oxygen- 
isotope effect on the pseudogap, in agreement with Eq.Sl 
Substituting SA = 16.5 meV and Srn^/mH = 4.3% into 
Eq. H we find W = 0.256 eV and t = 0.128 eV. Using 
m*^ = 2Ti'^ /to?, we calculate m*^ = 8.1 me, which is very 
close to that (8.3 me) inferred for the slightly underdoped 
YBaaCuaOe.BS with = 87.9 K (see ref. jsl]). 

From the A, Te, and W values of the ^^O and ^^O 
samples, we can directly calculate np{Tc) to be 0.0117 
and 0.0038 for the ^^O and samples, respectively. 
Substituting STc/Tc = -0.63% and SmU/mH = 4.3% 
into Eq. [71 we obtain dnp{Ta)/[x - np{Ta)] = -3.67%. 
With the np{Tc) values for the ^^O and samples, 
we find X = 0.227. Since the in-plane doping level x 
in the optimally-doped YBa2Cu3 06.95 was found to be 
0.264 (ref. our inferred doping level of 0.223 for 

HoBa2Cu408 is consistent with the fact that this com- 
pound is slightly underdoped. With m*l = 8.1 me and 
X = 0.227, we calculate XatiO) = 2248 A, which is close 
to Aa(0) = 2000 A for YBa2Cu408 (ref. [ll). 

The exponent a^. of the oxygen-isotope effect on 
m*^ is calculated to be 0.46 for the slightly underdoped 
YBa2Cu307_y and HoBa2Cu408. Since the exponent 
a*-^" of the copper-isotope effect on Tc is similar to a'~' 
in the underdoped regime 37|, |38| , we expect that a^. 
~ a^* and am* — 2a^* • The mass enhancement factor 
is then equal to exp(4a^.) = 6.3. The bare bandwidth 
D = 2W^exp(4a^j.) = 3.2 eV, in quantitative agreement 
with the band-structure calculations (39| . 

In the absence of charge localization, as the case of the 
stoichiometric HoBa2Cu408, the Bose-Einstein conden- 
sation temperature Tc is given by [Til] 



knTr 



2tnb{Tc) 



1 + \n{kBTc/2tc) ' 



(9) 
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where tc is related to the out-of-plane bipolaron mass 
m** as tc = /2m**cP. The out-of-plane bipolaron mass 
m** is deduced to be 518 rrie using rn** /rn**^^ = 64 at 
(ref. [3^]) and m*^ = 8.1 TOg- Then we calculate tc = 
0.160 meV. Substituting t = 0.128 eV, tc = 0.160 meV, 
and nh(Tc) = 0.1077 into the above equation, we cal- 
culate Tc — 78.9 K, in quantitative agreement with the 
measured Tc (79.0 K). We further show that (see Supple- 
mentary Information ) Eq. [201 can also quantitatively 
explain the underdopcd Lai.goSro.ioCu04 with Tc = 29 K 
but overestimates the Tc's of some optimally doped sam- 
ples. This implies that optimally doped cuprates are in 
the crossover regime from bipolaronic real-space pairing 
to the Cooper pairing of polarons [2^ . 

Apart from the striking isotope effects explained quan- 
titatively here there is abundant independent evidence in 
favor of bipolarons and the Bose- Einstein condensation in 
underdoped cuprate superconductors. In particular, the 
parameter-free estimates of the Fermi energy using the 
magnetic-field penetration depth 41| and the magnetic 
quantum oscillations yielded a very low value (below 
50 meV) supporting the real-space pairing in underdoped 
cuprate superconductors. Magnetotransport and ther- 
mal magnetotransport data strongly support preformed 
bosons in cuprates. In particular, many high-magnetic- 
field studies revealed the non-BCS upward curvature of 
the upper critical field Hc2(T) as a function of tempera- 
ture [43| as predicted for the Bose-Einstein condensation 
of charged bosons in the magnetic field [31 . The Lorenz 
number differs significantly from the conventional Som- 
merfeld value of the standard Fermi-liquid theory because 
the carriers are double charged bosons [i^. Direct mea- 
surements of the Lorenz number using the thermal Hall 
effect just above Tc jSl produce its value, which is about 
the same as that predicted by the bipolaron model. The 
unusual normal-state diamagnetism uncovered by torque 
magnetometery has been convincingly explained as the 
normal state (Landau) diamagnetism of charged bosons 
[131 ■ Single polarons, localised within an impurity band- 
tail, coexist with bipolarons in the charge-transfer doped 
Mott-Hubbard insulator, where the chemical potential is 
pinned within the charge-transfer gap due to bipolaron 
formation. This band-tail model accounts for two energy 
scales in ARPES and in the extrinsic and intrinsic tun- 
nelling, their temperature and doping dependence, and 
for the asymmetry and inhomogeneity of extrinsic tun- 
nelling spectra of cuprates jisf . 
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Supplementary Information— Isotope effects in 
high-Tc cuprate superconductors: Ultimate proof 
for bipolaron theory of superconductivity 

1 Microscopic derivation of the low-energy band- 
structure 

In highly polarizable ionic lattices like cuprate super- 
conductors both the Coulomb repulsion and the Frohlich 
electron-phonon interaction (EPI) are quite strong (of 
the order of 1 eV) compared with the low Fermi energy 
of doped carriers because of a poor screening by non- or 
near-adiabatic carriers In those conditions the BCS- 
Eliashberg theory P] breaks down because of the pola- 
ronic collapse of the electron bandwidth 3] so that one 
has to apply a non-adiabatic small polaron theory 

Here we sketch the microscopic derivation of the low- 
energy band structure. Fig. 1, using the analytical multi- 
polaron theory in the strong-coupling regime for highly 
polarisable lattices (more details are found in Refs. [4|). 

Quantitative calculations of the interaction matrix el- 
ements can be performed from pseudopotentials using 
the density functional theory (DFT) [H^. On the other 
hand, one can express the bare Coulomb repulsion and 
EPI through material parameters rather than computing 
them from first principles in many physically important 
cases 0. In particular, for a polar coupling to longitu- 
dinal optical phonons (the Frohlich EPI), which is the 
major EPI in polar crystals, both the momentum depen- 



6 



dence of the matrix element, M(q), and its magnitude 
are well known, |Af (q)| = 'y{q))huJo/\/2N with a dimen- 
sionless 7(g) = ^J'^Ke^^jK^)^ujjQ(p ^ where f2 is a unit cell 
volume, N is the number of unit cells in a crystal, is 
the optical phonon frequency, and k = eooeo/(eo — Coo)- 
The high-frequency, and the static, ep dielectric con- 
stants are both measurable in a parent polar insulator. 

The dielectric response function of strongly correlated 
electrons is apriori unknown. Hence one has to start with 
a generic Hamiltonian including unscreened Coulomb 
and Frohlich interactions operating on the same scale 
since any ad-hoc assumption on their range and relative 
magnitude might fail. 



2 eoolm - n| 



huJofii [u{m, q)dq + H.c] + Hj 



ph • 



(10) 



are essential in the expansion of the hopping operator 
aij. To deal with this challenging problem let us single 
out the coherent hopping in Eq. ljlip averaging aij with 
respect to the phonon vacuum, and consider the remain- 
ing terms as perturbation, H = Hq + Hp^ph- Here 

Ho = - ^^{UjSss' + (j-Sij)clcj + Hph (13) 

describes free phonons and polarons coherently propagat- 
ing in a narrow band with the exponentially diminished 
hopping integral, tij = T(m — n) exp[—g^{ni — n)], 

g^(m) = -^^7(g)'[l-cos(q-m)], (14) 



and 



H, 



p—ph 



I] )'Jss'c\cj 



(15) 



Here Tij = r(m — n) is the bare hopping integral, jj, 
is the chemical potential, i = m, s and j = n, s' in- 
clude both site (m, n) and spin (s, s') states, M(m, q) = 
(2A^)^^/^7((j) exp(iq • m), q, dq are electron and phonon 
operators, respectively, hi = c\ci is a site occupation 
operator, and Hph = X]q ^'^o('^q'^q + 1/2) is the polar 
vibration energy. 

In highly polarisable lattices with eo ^ 00 
the familiar Lang-Firsov (LP) Q canonical trans- 
formation e'^ is particulary instrumental with S = 
— Y^^^hilu^ai.q^dq — H.c.]. It shifts the ions to new 
equilibrium positions changing the phonon vacuum, and 
removes most of both interactions from the transformed 
Hamiltonian, H = e^He~^ , 



H 



i,3 



[cFijSss' + fl'5ij)c\cj + Hph, 



(11) 



where ct^ = r(m — n)xlXj is the renormalised hopping 
integral involving the multi-phonon transitions described 
with Xi = exp J2q u(ni, q)dq — H.c. , and jl = ji + Epis 
the chemical potential shifted by the polaron level shift. 



Ep — 



2ne' 



BZ 



(27r)3g2 



(12) 



Here, the integration goes over the Brillouin zone (BZ) 
and Ep = 0.647 eV in La2Cu04 [l[. The electron-phonon 
coupling constant is defined as A = 2EpN{0). In the 
case of 2D carriers with a constant bare density of states, 
iV(0) = nia^ /2iTh^ per spin, Ea. ([T^ places cuprates in 
the intermediate to strong-coupling regime, A ^ 0.5, if 
the bare band mass m > me (here a is the in-plane lattice 
constant). 

The number of virtual phonons in the polaron cloud is 
large in oxides and some other polar lattices, Ep/Tiujo > 
1 with the characteristic (oxygen) optical phonon fre- 
quency hujo ^ 80 meV, so that multi-phonon vertexes 



is the residual polaron-multiphonon interaction, which is 
a perturbation at large A. In the diagrammatic tech- 
nique the corresponding vertexes have any number of 
phonon lines. The second-order in Hp^ph polaron self- 
energy (Ep « —Ep/2z}?) and the phonon self-energy 
(Ep/i « —xTiLOo/zX^) are small, if A » l/\/2z ^8|](here 
z is the lattice coordination number and x is the atomic 
density of carriers). Hence the perturbation expansion 
in 1/A is applied. Importantly there is no structural in- 
stability in the strong coupling regime since |Sp;i| <C huo 

i- 

The LF transformation, Eq. (|TT|) is exact for any adi- 
abatic ratio huo/T{a). However, if the perturbation ex- 
pansion in 1/A is restricted by lowest orders, then it sig- 
nificantly overestimates polaron masses in the adiabatic 
regime, hujo/T{a) < 1, for the case of the short-range 
(Holstein) EPI (here T(a) is the nearest- neighbor bare 
hopping integral). The polaronic band narrowing factor, 
exp(— 5^) becomes very small for this EPI in the strong- 
coupling regime, which would eliminate any possibility of 
high temperature superconductivity and even metallicity 
of the small Hosltein polarons. 

However in the case of the long-range (Frohlich) EPI, 
Quantum Monte-Carlo simulations [9| show that the 
LF transformation provides numerically accurate polaron 
masses already in the zero order of the inverse- coupling 
expansion both in the adiabatic regime as well as in the 
non-adiabatic one for any strength of the Frohlich EPI. 
Moreover, such small polarons [9| and small bipolarons 
[lo| are perfectly mobile in the relevant range of the cou- 
pling and the adiabatic ratio. 

The perturbation Hp^ph has no diagonal matrix el- 
ements with respect to phonon occupation numbers. 
Hence it can be removed from the Hamiltonian in 
the first order using a second canonical transforma- 



tion n = e'^'^He-^'^ with (5*2 



Uj)c\cj\n) I {Er,> 



En), 



where En , En 



and In) 



are 
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the energy levels and the eigenstates of Hq, respectively. 
Taking into account that the polaron Fermi energy is 
small compared with the phonon energy at strong cou- 
pling and/or sufficiently low doping [ij, one can neglect 
the polaron contribution to — En ~ hujQ n'^ and 
project the second-order in 1/A Hamiltonian H onto the 
phonon vacuum |0) with the following result 

n = -"^{U^Sss' + ilS,j)clcj- 

^ ] cJ^sCnsCjjj/^/Cn's' , (16) 



mnm n ,ss 



where 



ym n _ -rp rp 



dte 



-St I 



[X}{t)x,{t)-1\XIX,,\Q), 
(17) 

and xl{t) is the Heisenberg multi-phonon operator ob- 
tained by replacing dq in X^ with exp(ia;ot). Calcu- 
lating the integral, Ea. (IT7|) with S +0 yields 



tijtj.y /(mn, m'n')'' 



Tiujo 



k=l 



klk 



(18) 



where /(mn, m'n') = (1/2A^) X;q7(9)^[cos(q-(m-n')) + 
cos(q • (n — m')) — cos(q • (m — m')) — cos(q • (n — n'))]. 

All matrix elements, Eq. ((T5)) . of the polaron-polaron 
interaction are small compared with the polaron kinetic 
energy except the exchange interaction, Jp{m — n) = 
such that /(mn, m'n') ~ 2g^(m — n). Using 
E^i2/Vfc!fc = -C - ln(y) + Ei*{y) with C w 0.577 
and Ei*{y) w e^/y (for large y) one obtains a substantial 
Jp(m) — T^(m)/25^(m)?icjo, which is larger than the 
nearest-neighbour polaron hopping integral, t{a)/Jp oc 
2huJog^e~^ /T{a) < 1. Keeping only this exchange we 
finally arrive with the polaronic "t-Jp" Hamiltonian 

hi 

+ ~ n) ('5m ■ Sn + ^n^nnj , (19) 



where = (1/2) E 



s,s' ^ms ' ss' ^ms' 



is the spin 1/2 op- 
erator (t are the Pauli matrices), JT-m = Tlis^i^ ^'^'^ 
fl = (1 + '^^Jp{i^) is the chemical potential further 
renormalized by Hp^ph- 

There is a striking difference between this polaronic 
t-Jp Hamiltonian and the familiar t-J model derived 
from the repulsive Hubbard U Hamiltonian in the limit 
U ^ t omitting the so-called three-site hoppings and 
EPI Til]. The latter model acts in a projected Hilbert 
space constrained to no double occupancy. On the con- 
trary in the polaronic t-Jp Hamiltonian, Eq. (I19|) there 
is no constraint on the double on-site occupancy since 
the Coulomb repulsion is negated by the Frohlich EPI. 
The polaronic hopping integral t{a) leads to the coher- 
ent (bi)polaron band and the antiferromagnetic exchange 



of purely phononic origin Jp bounds polarons into small 
superlight inter-site bipolarons. Last but not least the 
difference is in the " sign in the last term of Eq. 
proportional to fimfim which protects the ground super- 
conducting state from the bipolaron clustering, in con- 
trast with the "-" sign in the similar term of the stan- 
dard t-J model, where the phase separation is expected 
at sufficiently large J 12 1. 

The polaronic t-Jp Hamiltonian, Eq. (jlOp is analyti- 
cally solvable in the limit of sufficiently low atomic den- 
sity of carriers Q. Neglecting the first term in H, which 
is the polaron kinetic energy proportional to t{a) < Jp, 
one can readily diagonalise the remaining spin-exchange 
part of the Hamiltonian. Its ground state is an ensem- 
ble of inter-site singlet bipolarons with the binding en- 
ergy Ah = Jp localised on nearest neighbor sites. Such 
small bipolarons repel each other and single polarons via 
a short-range repulsion of about Jp. 

The kinetic energy operator in Eq. (jl9|) connects sin- 
glet configurations in the first and higher orders with 
respect to the polaronic hopping integrals. Taking into 
account only the lowest-energy degenerate singlet config- 
urations and discarding all other configurations one can 
project the t-Jp Hamiltonian onto the inter-site bipola- 
ronic Hamiltonian using the bipolaron annihilation oper- 
ators = 2"l/2(CintCm+ai - CmiCm+af), whcrC a COU- 

nects nearest neighbors jlO] . Such inter-site bipolarons 
are perfectly mobile since they tunnel via single-polaron 
transitions [lol . [isj . At finite temperatures single po- 
larons, thermally excited above the pseudogap, coexists 
with these bipolarons as shown in Fig. 1 of the main text. 

Small bipolarons are hard-core bosons with the short- 
range repulsion and a huge anisotropy of their effective 
mass since their inter-plane hopping is possible only in 
the second order of the polaron hopping integral [14| . 
The occurrence of superconductivity in bipolaronic sys- 
tems is not controlled by the pairing strength, but by 
the phase coherence among the electron pair s below the 
Bose-Einstein condensation temperature [15|. 
2 Parameter-free calculations of the Bose- 
Einstein consensation temperatures of some 
cuprate superconductors 

In the absence of charge localization, the Bose-Einstein 
condensation temperature Tc is given by [3] 



(20) 



where t is the bipolaron half-bandwidth and tc is re- 
lated to the out-of-plane bipolaron mass m** as tc = 
fi^ / 2m*c* d'^ [d is the inter-plane distance). The above 
equation can be written in terms of the measurable pa- 
rameters such as the inplane penetration depth and the 



supercarrier mass anisotropy constant 7 



ksTc 



167r3e2A2,(0) 
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32Ae^MAg,(0),^. (21) 

It is worth noting that the Tc value calculated from 
Eq. [21] should be slightly overestimated due to the fact 
that nb{Tc) is slightly lower than nf,(0) = x/2 . For 
Lai.9oSrpoCu04+a, Aab(O) ^ 291 nm (ref. fBl), 7 = 
43 (ref. [l8[), and x — 0.1. These parameters lead to Tc 
= 31.8 K, in quantitative agreement with the measured 
value of 29 K. For the optimally doped YBa2Cu306.95 
with T; = 93 K, Aafc(O) = 1600 A (ref. [11), x = 0.264 
(ref. H^l), and 7 = 8 (ref. (llj), so is calculated to be 
162 K, which is higher than the measured value of 93 K 
pointing to the BEC-BCS crossover Q. 
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